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Abstract 

Asymmetric nuclear matter is studied in the frame of relativistic mean-field theory, using scalar- 
isoscalar a, vector-isoscalar uj meson together with their selfinter actions, vector-isovector p meson 
with its cross-interaction with uj meson too, and scalar-isovector 5 meson as degrees of freedom. The 
model is used to parameterize the nuclear matter properties results calculated by more fundamental 
Dirac-Brueckner-Hartree-Fock theory and thus to provide an effective DBHF model applicable also 
to finite nuclei. Vector uj-p cross-interaction seems to be an useful degree of freedom for describing 
of the asymmetric nuclear matter, mostly due to its impact on density dependence of the symmetry 
energy. 
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I. INTRODUCTION 



The study of nuclear matter - hypothetical uniform infinite system of nucleons interacting 
via strong forces - has already been an important part of the nuclear physics development 
for several decades. It is a good starting point for both, nuclear physics of finite systems (e.g. 
structure and properties of finite nuclei, dynamics of heavy-ion collisions) and astrophysics 
(e.g. structure and evolution of stars). The early attempts were based on the nonrelativistic 
Brueckner-Hartree-Fock (BHF) theory (see e.g. Refs. |l| for review and references therein), 
and related the bare nucleon-nucleon (NN) interaction to nuclear ground-state properties 
in a parameter-free way with a limited success. 

The break-through was achieved when the relativistic extension of the BHF theory (so 
called Dirac-Brueckner-Hartree-Fock (DBHF) approach) was developed ODD and suc- 
cessfully applied to nuclear matter problems. An essential feature of the DBHF is the in- 
corporation of the relativistic dynamics, governed by the Dirac equation with strong scalar 
and vector fields. The explicit treatment of the lower components of Dirac spinors gives a 
rise to strongly density-dependent relativistic effects. They shift the nuclear matter satura- 
tion points ("Coester band") towards empirical values Subsequently, a great effort has 
been devoted to the successful DBHF description of both, symmetric and (to lesser extent) 
isospin-asymmetric nuclear matter. 

Thus, the DBHF approach is currently considered to be a microscopic parameter-free 
nuclear model based on realistic NN interaction. However, due to its complexity, this so- 
phisticated approach is successfully manageable for nuclear matter properties only; finite 
nuclei are at present beyond the scope of this model. To overcome this restriction, several 
approaches were developed which relate the DBHF output for nuclear matter to the param- 
eters of the relativistic mean-field (RMF) theory ((JO]. The RMF approach is widely used 
and powerful phenomenological tool for various aspects of nuclear many-body problems 
which provides an effective framework for calculation of both, nuclear matter and (contrary 
to DBHF) also finite nuclei. The nonlinear RMF approach 8j has already been proven to be 
a reliable tool for the calculation of normal nuclei close to the valley of stability jjj, exotic 
and superheavy nuclei |ic| . 

Now new experimental facilities are available to study properties of exotic nuclei with high 
isospin asymmetry. Additionally, increasingly more precise observations and measurements 
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of properties of neutron stars and supernovae have been carried out. This naturally brings 
a need for better description of isospin degree of freedom, which can be done by enhancing 
;he isovector meson sector. The isovector scalar 5 meson and vector cross-interactions 



12j were included into RMF for this purpose. 

The goal of this paper is obtaining of effective parametrization applied to asymmetric 
nuclear matter, using different degrees of freedom in order to study influence of 5 meson 
and vector meson cross-interaction on quality of the reproduction of DBHF results, as well 
as their influence on calculated nuclear matter properties, especially on density dependence 
of nuclear symmetry energy. 



II. THEORETICAL FRAMEWORK 



The starting point of the model is Lagrangian density that introduces nucleon field 
isoscalar-scalar meson field a, isoscalar- vector meson field w, isovector- vector meson field p 
and isovector-scalar meson field 5 (pion field does not contribute, because it is pseudoscalar, 
and nuclear matter is parity invariant), and takes a form 



£(V>, <r, uj, p,5)=i/j bv(^ - g^ - (M - g a a)\ if; 

-\b a M{g a af - \c a {g a af + \c^gl^f 
+ \{d^5 - m}5 2 ) + \vP? p p^ - \p^.p^ 
+\k v {glp^)(glu^) 

-gpp^i^T"^ + gs54>Tip , (l) 

where antisymmetric field tensors are given by 



P^v = d u p^- d^p v , 

and the symbols used have their usual meaning. The parameters entering the lagrangian 
are M that denotes the nucleon free mass, whereas m a , m w , m g and m$ are masses assigned 
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to the meson fields. The first term together with the last two ones describe interaction 
of isoscalar and isovector mesons with nucleons where the strength of these interactions is 
determined by dimensionless coupling constants g a , g w , qv, and gs- Three terms in the third 
ine rep resent cubic and quartic scalar selfinteractions [8| and quartic vector selfcouplings 
kJQI) t ne strength of which is also given by dimensionless selfinteraction coupling constants 
b a , c a and c w . The second and fourth lines represent free (noninteracting) lagrangian for all 
mesons, and the fifth line realizes cross-interaction between u and p mesons characterized 
by cross-coupling constant Ay 

Constraint of stationarity of the action leads to the well-known Euler-Lagrange field 
equations and equations of motion follow after their application to the lagrangian (JTJ. This 
produces the Dirac equation for nucleon field 

[7„(i0" - g^ - g pPp .r) — (M — g a a - g s S.r)] i> = . (2) 

Isoscalar meson field a, to are then described by Klein-Gordon and Proca equations, 
respectively, 

(d^ + ml)a = g a $i> - b a M(g a a) 2 - c a {g a af] , (3) 

-Avfiprffguutp] . (4) 
Analogically, isovector p and 5 meson fields read, 



d^pT + my = g p $YTi(> - k Y g p p,glu>^\ , (5) 

{d^ + m 2 s )5 = g s ?pT?P . (6) 

Due to the fact that these equations are nonlinear, nowadays there is known no suitable 
method to solve them. The way to avoid this is to replace the operators of meson fields by 
their expectation values - the so called mean-field approximation. The fields are thus treated 
as classical c-numbers. Its reasonability increases with increasing baryon density. The second 
approximation introduced is the non-sea approximation which doesn't take account of the 
Dirac sea of negative energy states. 
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In this model we are dealing with static, homogenous, infinite nuclear matter that al- 
lows us to consider some other simplifications due to translational invariance and rotational 
symmetry of nuclear matter. This causes the expectation values of space-like components 
of vector fields vanish and only zero components - po and ujq - survive |9]. In addition, ro- 
tational invariance around third axe of isospin space results in taking into account only the 
third component of isovector fields - p^ and 5^ jl, 11]. The above-mentioned can formally 
be written as 



a — > (a) = a , 

Ptx ► (P/i) = Po } , 

5 — > (S)= 5 (3) . 

Having inserted the above simplifications the field equations are reduced and we can 
easily obtain potentials of both isoscalar meson fields 



U a = -g a a = -&\p s ~ KM(9^) 2 

-Ca(g*v) 3 } , (7) 
g 2 s 

U w = QujUq = -^[p B - c u (g u u ) 
mi 

-U 2 p A v (g w iu )] , (8) 

and isovector meson fields 



Up ee g p p { Q ] = -^[ip7°T 3 ip - gpP^Avig^o) 2 } = 
m p 

%4-l) PB -g p p^A v U 2 }, (9) 
m p j\ 

Us == -gsS {3) = -4^ = 4(Pn - P S P ) , (10) 
where scalar density ps is expressed as the sum of proton (p) and neutron (n) part 

PS = (H>) =Pp+Pni (11) 
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which are given by 



2 r k i M* 
of = o/ d 3 k - T r,i = p,n. (12) 

(27T) 3 i0 (fc2 + Mf) 1/2 

In the equation (fT2"j) hi is nucleons' Fermi momentum and M* M* denotes proton and 
neutron effective masses, respectively, which can be written as 



M* p = M - g a a - g 5 m , (13) 
M: = M - g a a + g 8 6^ . (14) 

One can see that condensated scalar a meson field generates a shift of nucleon mass, in 
consequence of which nuclear matter is described as a system of pseudonucleons with masses 
M* moving in classical fields a, loq and /?o j where additionally 5 meson field is responsible 
for splitting of proton and neutron effective masses, which is an important feature of 5 
meson influence on the nuclear matter saturation mechanism and its properties. The S 
meson seemed to be an useful degree of freedom in describing of asymmetric nuclear matter, 
indicated by its influence on e.g. stiffness of equation of state, slope and curvature of 
symmetry energy and properties of warm asymmetric nuclear matter |lll . Ilq . 

The solution requires to be performed selfconsistently, which can be clearly seen from 
equations (J7j)- (fTU)l where a potential (J7|) must be solved using iterations. 

The baryon density is given by 

ffi =K*>^f (15) 

with kp beeing an average Fermi momentum. It can be seen that scalar density (jllj) is less 

1 /2 

than baryon density due to term M* / (k 2 + M* 2 ) that causes reduction of the contribution 
of rapidly moving nucleons to scalar source term. This mechanism is responsible for nuclear 
matter saturation in the mean field theory and essentially distinguishes relativistic models 
from nonrelativistic ones. 

Cross-coupling of the u and p mesons requires also selfconsistent calculation of equations 
(jHJ) and 0, with iterative procedure for u potential. 

By reason that 5 field splits nucleon effective masses, the proton and neutron Fermi 
momenta will be also splitted, while they have to fulfil 
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pB = Pp + pn = ——3 dk + —— s dk, 16 
(27T) jo (27r) Jo 

where &p is average Fermi moment of the matter, k p , k n are Fermi momenta of protons 
and neutrons, respectively. The different values of Fermi momenta have consequences for 
transport properties of asymmetric nuclear matter. 

To obtain formula for energy density of nuclear matter it is essential to have cognizance 
of the energy tensor, in continuum mechanics defined as 

d$j dC 

- -a^ + dxU d ( d $./ dXii ) > y 17 > 

where $j generally denotes physical fields. The energy density of such a system is the zero 
component of the energy tensor e = (Too), and finally the binding energy per nucleon is 
related to energy density by 

E b = — -M. (18) 
Pb 

The energy density per nucleon is a starting quantity for further properties of nuclear 
matter. Incompressibility is given as its second derivative with respect to baryon density, at 
the saturation point 



K = 9 



2 d\e 



P 



(19) 

J P=PO 



dp 2 p B 

Symmetry energy of nuclear matter is defined as a second derivative of binding energy 
per nucleon with respect to the asymmetry parameter a = (p p — p n )/(p p + p n ): 

e{p, a) = e(p, 0) + S 2 {p)a 2 + S A (p)a\ (20) 
where parameters S 2 , S4 are defined as 



S 2 

S 4 ~- 



1 

' 2 
1 

24 



d 2 e(p, a) 

da 2 
~ d 4 e(p,a] 

da 4 



a=0 



a=0 



(21) 
(22) 



Parameter S 2 is often used as symmetry energy itself, argumented by negligible contri- 
bution of higher order S4 parameter, especially for densities relevant for common nuclei. 



III. RESULTS AND DISCUSSION 



The mean-field parametrizations were obtained by calculation using three different DBHF 
results for nuclear matter as initial data - results of Li, Brockmann and Machleidt jl?! ]. 
results of Lee, Kuo, Li and Brown IllJ an d finally calculations of Huber, Weber and 



Weigel 



20 



2l|. 



Based on the realistic and relativistic NN interaction of the Bonn group, in the work 



171 ] there were performed DBHF calculations that yield an effective NN interaction, and 
subsequently the single-particle potentials, equations of state, nucleon effective masses, and 
speed of sound for both symmetric and neutron matter were studied. The Bonn A potential 
reproduced quantitatively the empirical saturation properties of nuclear matter as well as 
nucleon effective mass. 

The work deals with asymmetric matter also using the DBHF approach with Bonn 
A one-boson-exchange NN interaction. Not only saturation properties, but in addition 
even the empirical value of the symmetry energy at the saturation density were reproduced 
satisfactorily. Isoscalar meson potentials for symmetric matter authors calculated in |ig| . 
inally energy per nucleon for several asymmetries using DBHF approach is calculated 



in j2pj, together with proton and neutron scalar and vector potentials from [2jJ making 



possible to fit also the mean-field parameter for coupling of 5 meson to nucleons. Further- 
more, in this two works it is tested also momentum dependence of self-energies, however, 
due to momentum independence of mean-field theory coupling constants, the momentum 
independent results were chosen to fit. 

As it is easily seen from equations for meson potentials and for energy per nucleon, the 
squares of coupling constants appear exclusively in ratios with meson masses and thus one 
can fix the meson masses to experimental values without any physical restriction of the 
RMF. The meson masses considered in this work are m a = 550 MeV, m u = 783 MeV, 
m p = 770 MeV and m s = 980 MeV. 

The first fit was performed using a, to, p mesons and scalar isoscalar cubic and quartic 
selfinteractions, vector isoscalar quartic selfinteractions as well as cross-interactions between 
vector mesons, fitting energies per nucleon for symmetric and neutron matter and symmetric 
isoscalar potentials. The corresponding parameter set obtained is listed in Table H] - for 
comparison with (parametrization MA) and without (MB) vector meson cross interaction. 
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TABLE I: Parameter sets resulting from the RMF fit to DBHF results of Machleidt et al. [1 





MA 


MB 


al 


106.85 


112.27 


al 


180.61 


204.36 


$ 


18.445 


9.4932 


K 


-0.0025823 


-0.0029820 


Ca 


0.011529 


0.013345 




0.015849 


0.020449 


Av 


0.25857 




X 2 /N 


2.76 


9.95 



This parametrization reproduces the DBHF results satisfactorily in whole fitting range of 
densities relevant for common nuclei, which can be said both for energy per nucleon and also 
for isoscalar a and to meson potentials. Growth of the p potential is decreasing with baryon 
density, being result of the isovector cross-interactions. This has impact on the energy per 
nucleon especially for extreme asymmetries and consequently on curvature of symmetry 
energy, which is plotted in the Fig^ upper panel. The cross-interactions significantly affect 
density dependence of the symmetry energy - they increase rise of symmetry energy below 
0.24 fm -3 and slower it above this density. Symmetry energy at the saturation density is 33.3 
MeV, which is in accordance with experimental value of about 34 MeV. The incompressibility 
of symmetric matter at the saturation density is 347 MeV. 

For better description of asymmetry behavior of matter it is of course profitable to use not 
only symmetric and neutron case, but also other nonsymmetric cases with partial fraction of 
protons, as was calculated with DBHF theory by Lee et al. 18, Q|. Results of the fit with 
the same degrees of freedom as were in the previous case, performed for binding energy per 
nucleon for several asymmetries as well as for symmetric matter scalar and vector potentials, 
are listed in Table ITT1 (parametrization LA), also with parameter set obtained using model 
without vector meson cross-interactions (LB). All of the relevant physical quantities for all 
of the asymmetries are reproduced closely. Also in this parametrizations there is strong 
influence of cross-interaction, which is displayed in Fig Q lower panel. There is increas- 
ing growth of symmetry energy below approximately 0.25 fm -3 and deceleration in higher 
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FIG. 1: Density dependence of the symmetry energy for two different parametrization sets. Upper 



Q. 



panel shows results for the fit of Machleidt et al. DBHF results [Tj| , with mean-field parametriza- 
tions listed in Table QJ where a, u mesons with their selfinteractions and p mesons was used as 
degrees of freedom, and both with and without inclusion of vector meson cross-interactions (VCI), 
plotted with solid and dotted lines, respectively. The lower panel displays analogical calculation 



results for the fit of Lee et al DBHF results 



with parametrizations in Table ITTI 



density region. This fact will be comment more closely in the last data set. 

In the theoretical framework we was dealing also with isovector-scalar sector of nucleon- 
nucleon effective interaction. In the work [20] there was calculated binding energy per 
nucleon for several asymmetries in the DBHF approach, using Bonn potential B with density 
independent self-energies. Additionally in 2l( authors performed calculation also for proton 
and neutron scalar and vector potentials, enabling us to fit not only p potential but also 
S potential value. The fit was thus using a, u mesons, their selfinteractions, p meson with 
its cross-interaction to u meson, and finally even 5 meson as degree of freedom. Results of 
the fit for binding energy per nucleon for several considered asymmetries are drawn in the 
Fig |21 Fit values (represented by lines) follow closely the DBHF results (scatter symbols) 



10 



TABLE II: Parameter sets resulting from the fit of Lee et al. DBHF calculations 18L llS |. 





LA 


LB 




103.91 


102.11 


gl 


147.84 


146.73 


$ 


17.432 


9.6697 


K 


0.00097186 


0.00083559 


Ca 


0.0012694 


0.0012411 




0.0054204 


0.0051878 


Av 


0.18790 




X 2 /N 


1.69 


2.62 



TABLE III: Parameter sets resulting from the fit of Huber et al. DBHF results 2C 



HA 



HB 



HC 



HD 



9l 

gl 



9S* 

K 

X 2 /N 



90.532 
108.95 
36.681 
28.739 
0.0043852 
-0.0052045 
-0.0001421 
0.10647 
2.05 



86.432 
106.89 
28.795 
25.170 
0.0033779 
-0.0037762 
-0.0010509 

3.80 



91.110 
109.26 
20.804 

0.0044388 
-0.0052076 
-0.0000385 
0.34805 
5.85 



87.591 
107.61 
15.335 

0.0035745 
-0.0039753 
-0.0007753 



for all asymmetries. Corresponding parameter set is listed in Table ILTTl fHAL together with 
fit parameters obtained without vector cross-interactions (HB). 

Negative values of the quartic selfinteraction constants are consequence of the fact, that 
dependence of symmetric scalar and vector potentials ensued from DBHF calculations on 
the baryon density is almost linear, which brings some difficulties into the determination 
of the selfinteraction force of isoscalar mesons, and thus results in negative and positive 
second derivative of vector and scalar potential, respectively. One of the explanations of 
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FIG. 2: Density dependence of the binding energy per nucleon for five different asymmetries from 
pure neutron matter to symmetric nuclear matter, resulting from mean-field theory fit (represented 



by lines) of Huber et ah DBHF results 



20 



2l| (scatter symbols), with corresponding parametriza- 
tions listed in Table ITTTI Isoscalar a and u> mesons with their selfinteraction, isovector p mesons with 
(solid lines, parametrization HA) and without (dotted lines, parametrization HB) cross-interaction 
with to mesons, and 5 mesons was used as degrees of freedom. 



this verity could be the fact, that the Van-Hove theorem |22j is, unlike for the DBHF theory, 
fully consistent with the mean-field model only. The fit of Lee et al. DBHF calculations 
is not the case (see Table HTl - all selfinteraction constants have positive sign), due to more 
distinct curvature of density dependence of the potentials (not shown in this paper). The 
incompressibility of nuclear matter is 235 MeV. 

Simultaneously with binding energies fit to the symmetric scalar and vector potentials 
(Fig |3J) and also to above-mentioned proton and neutron scalar and vector potentials for 
proton fraction 0.125 (Fig H} was performed. As can be seen all of these potentials are 
reproduced satisfactorily within several MeV. 

Evaluation of p-u cross-interaction influence on symmetry energy and its comparison 
with 5 meson influence is possible from Fig There it is drawn density dependence of 
the symmetry energy, where each of the lines is correspondent to different degrees of free- 
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potentials of symmetric nuclear matter, resulting from the fit of Huber et al. 
, with the same parametrization and notation as in the Fig [21 




FIG. 4: Scalar and vector potentials of protons and neutrons of the asymmetric matter with proton 
fraction Z/A = 0.125, for the parametrization and notation identical with previous figures. 
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FIG. 5: Density dependence of the symmetry energy of nuclear matter, as resulted from the fit of 
Huber et al. DBHF results I20I . 21 1 for several considered degrees of freedom (all models include 
also isoscalar scalar and vector selfinteractions, see Table 1111}) : HA: a,uj,p and 5 mesons and p-uo 
cross-interaction; HB: cr,uj,p and 5 mesons; HC: cr,LV,p mesons and p-to cross-interaction; HD: a,LO 
and p mesons. The inset picture shows the symmetry energy behavior in a wider density region. 



dom used (see Table UTT]) : HA contains both cross- interactions and 5 meson contribution; 
then inclusion of the 5 meson without vector cross- interaction is accomplished by HB; HC 
considers the vector cross-interaction but without 5 meson, and finally basic model with a 
and uj mesons with their selfinteractions and p mesons only (HD). Vector cross-interaction 
has significant impact on symmetry energy, where its influence has analogical nature as 
in the previous two data sets, thus supporting those results also in the case of S meson 
inclusion. It increases growth of the symmetry energy in this case below approximately 
0.13 fm -3 , which is important for description of properties of exotic nuclei near the dripline 
and is similar as consequences of density dependence of the isovector couplings due to Fock 
contributions, which was calculated in j^. Above 0.13 fm~ 3 it slowers symmetry energy 
rise thus implicating an impact on higher density behavior of matter (high energy beam 
collisions, supernovae explosions, neutron stars properties), and it is also in agreement with 
recent Hartree calculations j^j], where similar high-density influence of cross-interactions 
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was concluded. In comparison 5 meson contributes relatively slightly to behavior of sym- 
metry energy - for lower and intermediate densities its contribution is almost inconspicuous 
thus not significantly affecting properties of atomic nuclei. For higher densities in absence 
of cross-interaction it slightly increases symmetry energy. Also this result is in concordance 
with recent calculations of other authors [25j. However, the presence of cross- interaction 
results beside strong curvature of symmetry energy also in contrary (compared to previous 
case) effect of 5 meson - for higher densities it softly fortifies decrease of symmetry energy 
growth. This leads to conclusion, that p-uj cross-interactions seems to be an important 
degree of freedom which should be used in further calculations. 



IV. SUMMARY 



In this work the relativistic mean field theory was used to obtain an effective parametriza- 
tion of the properties of asymmetric nuclear matter calculated by more fundamental Dirac- 
Brueckner-Hartree-Fock theory. The energy per nucleon together with the symmetric 
isoscalar potentials were fitted, and simultaneously also proton and neutron scalar and vec- 
tor potentials. Isoscalar a, uj mesons with their selfinteractions, and isovector p, 5 mesons 
with p-uo cross- interaction were used as degrees of freedom and parameters of the fit. Gen- 
erally a good reproduction of both the energy and the potentials was reached and thus the 
parameter sets are representing an effective DBHF description of asymmetric nuclear matter 
at normal baryon densities applicable as well for calculation of finite nuclei properties. The 
cross-interaction between p and u mesons turns out to improve reproduction of properties 
of asymmetric nuclear matter and additionally it increases symmetry energy in common nu- 
clei density region, e.g. approximately bellow saturation density, and decreases this energy 
above saturation point, thus having consequences for properties of finite nuclei, especially 
with large isospin asymmetry, and also for description of nuclear matter at higher densities, 
relevant in high energy nuclear collisions and several astrophysical processes and phenomena 
(e.g. neutron star properties and supernovae explosions). Isovector 6 meson also improves 
quality of the mean-field model, but without such a strong impact on density dependence 
of symmetry energy. This results imply that p — u cross-interaction is very useful for better 
description of nuclear matter and that it could be important for calculations of properties 
of finite nuclei with high isospin asymmetry. 
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